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1. Introduction 

Within perturbative QCD (pQCD) at high energies it is by now well understood that an ad- 
equate treatment should account explicitly for the transverse-momentum degrees of freedom of 
the participating partons. Specifically in the high-energy limit, this is accomplished by the so- 
called -factorization which ultimately goes back to works on the Balitsky-Fadin-Kuraev-Lipatov 
(BFKL) equation for the small-;c evolution of the unintegrated glue ([|I|], for recent reviews on k±^- 
factorization and more references see ^]). An underlying assumption of the phenomenology 
based on (linear-) ^^-factorization is that it is sufficient to consider only single gluon (reggeon) 
exchange in the parton-target interaction. Certainly, this can only be warranted in a finite energy 
range, although perhaps covering most of the small-jc domain oi Y*p collisions at HERA. Indeed, 
in a color-dipole picture [Q], HERA data on diffractive deep-inelastic scattering imply, that for 
color dipoles Oei/Ctot *C 0.5, that is unitarity constraints must be weak. Still, eventually the s- 
channel partial waves of single-gluon exchange will overshoot unitarity bounds, and the inclusion 
of multiple gluon/reggeon exchanges associated with absorptive/unitarity corrections is mandatory. 
Although some semi-quantitative phenomenology exists (e.g.[|5|]), a systematic theory for proton 
targets is absent. Here heavy nuclei as targets offer a welcome testing ground for unitarity effects. 
Multiple gluon exchanges are enhanced by the large thickness of a heavy target nucleus, and one 
can concentrate on those unitarity corrections which grow with the size of the target A. From a 
different viewpoint [^, the opacity of heavy nuclei generates a new large scale Q\ oc 7^(b), the 
saturation scale, where TAih) is the nuclear thickness (see e.g. for a review and references). 
Especially in view of present and future programs on forward physics at RHIC and LHC, the fate 
of kj^ factorization in a regime of strong nuclear absorption is an important issue. 

For illustration, consider the production of dijets in the virtual photon fragmentation region 
in 7*-proton collisions. Here the dominant contribution comes from the qq final states and can 
be viewed as a hard photon-gluon fusion Y'gt — > qq. While in the familiar collinear factorization 
the final state is a pair of back-to-back jets, in the lowest order of -factorization the quark 
and antiquark jets are decorrelated in azimuth, and the distribution in the decorrelation momentum 
maps out the proton's unintegrated glue 



(1.1) 



Here p is the transverse momentum of the quark-jet, z is the fraction of the photon's lightcone mo- 
mentum carried by the quark, and the transverse momentum K transferred by the gluon quantifies 
the azimuthal decorrelation of the <7^-system. Finally p) is the momentum space Ughtcone- 
wavefunction (LCWF) for the transition 7* — > qq, and f{x,Vi) K~'*5G(;c,K)/5log(K^) denotes 
the unintegrated gluon distribution of the proton target, which e.g. is directly proportional to 
the amplitude of hard diffractive dijet production [Q]. When addressing the same observable in 
deep inelastic scattering (DIS) off nuclear targets [[l^], using a nuclear unintegrated glue that ful- 



fills the same proportionality to hard diffractive amplitudes Jllp, we found that the linear 
factorization is completely broken. The spectrum of dijets off nuclei is an involved nonlinear 
functional with little resemblance^ of the free-nucleon result eq. ( |1.1| ). It incorporates distinct 

' See eq.(^j]) below. 
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physical effects, such as the coherent nuclear distortion of the qq LC-WF as well as incoherent 
transverse-momentum broadening through intranuclear multiple scattering, which arise from the 
same unitarity-driven dynamics. We dubbed the emerging momentum-space formalism Nonlinear 
Factorization, and subsequently developed the description of other relevant hard subprocesses. 



such as TT —> qq, q — > qg, g —* gg tg ^ qq [|12|, 13, 14]. In this contribution we will address the 



most recent developments: first, the evaluation of virtual radiative corrections to single-jet ob- 
servables, and the small-;c evolution of the nuclear unintegrated glue from the i^-channel unitarity 



relations [15], and second an extension to more refined observables, such as topological cross 
sections and the related unitarity rules [ p^ of Abramovsky-Gribov-Kancheli (AGK [p^) type. 
This will include prolegomena to a Reggeon Field Theory (RFT) interpretation of the nonlinear 
A: -factorization results. The breaking of linear /cj^ -factorization has also been reported in other 



approaches, we mention fllSQ. 



2. Nuclear unintegrated glue []I^ and its evolution [ jl^ 

Swiftly moving partons propagate along fixed impact parameter trajectories, and their impact 
parameter is conserved in the interaction with the target A. Consider a parton a in an arbitrary color 
representation, and let Sa(ba) be the S'-matrix of the (bare) parton-target interaction in the target 
rest frame at an impact parameter b^. Then, the S'-matrix for the ad color-dipole of size r = — b^ 
scattering at b 

c; .Hr^ (A|Tr[S,(b,)S^(b,-r)]|A) 

can be evaluated by standard Glauber-Gribov[[l^, ^ ] techniques at a boundary value a;^ = (2/?a'«a') ^ 
{Ra is the nuclear radius, mn the nucleon mass) as Sfla(b,r) = exp[— ^aaa(xA,r)7A(b)], and serves 
to define the nuclear unintegrated gluon distribution 0(b,;cA,p) 

^{h,XA,v) = J ^^S«<5(b,r)exp[-/pr] =exp[-VA(b)]5P)(p) + 0(b,XA,p). (2.2) 

Using its Glauber-Gribov relation to the free-nucleon color-dipole cross section a{xA,r)aa one 
derives a convenient representation of the coherent nuclear glue (j){h,XA,p) in terms of multiple 
convolutions /^■'^(p) = {f ^ f^^^^^)ip) of the free-nucleon unintegrated glue /(p): 0(b,:tA,p) = 
^wy(V/i (b) (p). Here the Poisson weights ^ wj = exp[— VA(b)]v{(b)/7! represent the prob- 
ability to find j overlapping nucleons at impact parameter b. The salient features of 0(b,A;A,p) 
can be briefly summarized: first, for realistic free-nucleon input and heavy nuclei there emerges a 
saturation scale Q^Ixa) ~ 0.8 1 GeV^, second, at large p^ we have a Cronin-type antishadowing 
enhancement that is calculable parameter-free, third it furnishes a linear /cx-factorization of inclu- 
sive DIS, forward single-jets in DIS and diffractive dijet amplitudes. To study evolution to small 



x-^xa, one has to follow the strategy of [[21|], and integrate out the effect of an additional s-channel 
gluon, schematically: 

5S„„-(x,b,r) oc log J \Yaag\^ {Saag " Sao) ■ (2.3) 



^The nuclear opacity VA(b) = ^CTo(^A)?A(b) involves the dipole cross section for large color dipoles Go{xa) and 
the nuclear thickness function 7a (b)- 
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VRC 



Figure 1: The real emission (REC) and virtual radiation (VRC) contributions to the spectrum of leading 
quark jets in inclusive production off nuclei. The (pseudo)rapidity of radiation correction gluons runs be- 
tween the quark jet rapidity and the nuclear boundary condition rapidity r\A = log(l /xa), where xa is defined 
in the text. 



This produces the first step of the small-;c evolution, a closed form of which does not exist, and is 
equivalent to the first step of the Balitsky-Kovchegov (BK) equation [22]. Using the relation 
one may translate (2.4) to momentum space. Concentrating on color-triplet partons (quarks), one 



would obtain [fi31] 



50(b,.x,p) =log( — ) (J^FA-L«'0)(b,XA,p) + ^[0](b,.XA,p) 



(2.4) 



Here we decomposed the nuclear evolution into the linear, BFKL-evolution piece, {J^bfkl ® 

uinely nonlinear component represented by the quadratic functional ^ 

^[^](b,;cA,p) = y"<i^qi/q2</'(b,.XA,qi){[^(p + qi,p + q2)-^(p,p + qi) 

-^(p,p + q2)]0(b,jCA,q2) -</'(b,JCA,p)[^:(qi,qi +q2 + p) -^(qi,qi +p)]}. 

The latter generates a higher-twist contribution to the large-p tail of the nuclear collective glue, 
but it does not exhaust the nuclear higher twists, recall that the boundary condition itself contains 
substantial - but an^Z-shadowing - higher twist contributions. 



3. Forward jets from the breakup of valence quarks [ )15[ J 

Consider now the production of jets/high px^ particles at forward rapidities, i.e. in the beam 
fragmentation region in a high energy proton-nucleus collision. The dominant beam partons are 
evidently valence quarks. The density of valence quarks is a steep function for ;c — > 1, and the pro- 
duction of high forward partons must be highly sensitive to their energy loss. Here we discuss 
the perturbative, radiative energy loss that originates from the transition q — > qg. The dynamics 
of partons with rapidities T] > tja = log(l/;cA), proceeds coherently over the whole nucleus (for 
some notation see Fig. |l|). Notice that in addition to the coherent excitation q — > qg, in the 

^here, for perturbative large p;, A^(pi,p2) = (pi — P2)^/PiP2- For ^ regularized form valid also at small p;, see 

0. 
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factorization one must also add the quasielastic scattering of the beam quark, where its transverse 
momentum is balanced by (mini-) jets with rapidities T] < T]^. To quantify the /^j^ -dependent ra- 
diative energy loss, one must equally evaluate the radiative correction to the quasielastic scattering, 
i.e. the virtual corrections in terms of fast r\ > r\A parton degrees of freedom depicted as VRC 
in Fig. The calculation of radiative corrections, real as well as virtual is most expediently per- 
formed in impact parameter space. To this end, one expands the LCWF of the beam quark into its 
bare quark and quark-gluon components. By the conservation of impact parameters, the S'-matrix 
acts on the physical quark state in a simple form: 

S\qpHys{h)) = [S,{h) + 5S,mqphys{h)) + [S,{h^^^^^ 

(3.1) 

Evaluation of two-parton and single-parton spectra from the excitation operator Sex = Sq^'q)Sg(!^g) ' 



^^(b) is discussed in detail in [[Tj, The important point is that Sex determines through the uni- 
tarity relation also the radiative correction SSqih) to the quasielastic scattering of the quark, which 
reads: 

5Sq{h) = J nb,,z,}\v,f[S,{Kj)SM-S,{h)] . (3.2) 

It is important to stress, that the LCWF/color-dipole S-matrix approach is fully general in that 
it applies equally to the single-gluon exchange encountered on the free nucleon as well as to the 
multigluon exchanges relevant to strongly absorbing nuclei. Also notice, that we do not require the 
multi-Regge type kinematics of Tjg ^ rjq, but allow for an arbitrary rapidity spacing of scattered 
quark and the radiative gluon as is appropriate for the evaluation of the radiative energy loss. We 
skip the detailed form of the single-parton spectra, as well as the case of incident color-octet 
partons, and close this section with a remark on the small-jc evolution properties of the single- 
quark spectra. At the boundary value xa one can easily check that the latter is directly proportional 
to the target unintegrated glue, dOQei{q q,XA)/d^hd^p = 0(b,;cA,p). The intriguing result is 
now that this relation holds, to the leading-log (LL(l/jc)) also after the radiative corrections, 
that is 

dOQei{q^q,x) \ , s^/u \ o^ 

^2b^2p = (j){h,XA,p) + 5(j){h,x,p), (3.3) 

with the same 5^{h,x, p) of eq.([2.4[) that defines the first LL(1 /x) iteration of the nonlinear small-x 



evolution of the nuclear unintegrated glue. We recall also the earlier result [13], that the spectrum 
of midrapidity gluons is linearly factorizable in terms of the nuclear unintegrated glue. All of this 
suggests a deeper role of the coherent nuclear glue in terms of a to be developed RFT description. 
Some of these aspects we shall address in the following. 



4. Topological cross sections and unitarity cutting rules fll^ ] 

In linear k± factorization the partonic transitions a ^ be of interest proceed via single gluon 
exchange and leave the recoiling target nucleon debris in a color octect state. There is an ob- 
vious connection between the unintegrated gluon distribution and a cut (BFKL-) Pomeron. For 
example, the quasi-elastic scattering of parton a with a color excited nucleon in the final state is 
dOQei/d^p = Ca/2Cf /(x,p), where Ca,CF are color Casimirs and the nuclear unintegrated glue / 
represents the coupling of the color-octet ^-channel exchange to the nucleon. The situation is more 
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involved for the heavy nucleus target. On the beam side, during its coherent interaction with the 
nucleus the be parton system will evolve over all possible color multiplets in the relevant product 



representation space. The at first sight formidable problem in color algebra can be solved [ |14| ] 
using a technique originally due to Zakharov, where by concentrating on the Z7cl?'c'-density matrix 
one deals with a four-parton system in an overall color-singlet state. The intranuclear evolution 
problem is then solved by a color-coupled channel four-body S'-matrix S{C) = exp[— ^r('^)(C)lA] 
(here C collectively denotes the relevant impact parameters). The building block is the color-dipole 
cross section operator r(^)(C) for the bcb'c' system for the free-nucleon target. It is a matrix in the 
space of possible color-singlet states \RR) = \ {bc)R ® (b'c')^). The overall-singlet property of the 
bcb'c'-system implies that its matrix elements are infrared safe combinations of the free-nucleon 
color dipole cross section. It is important to distinguish the following aspects of color channel cou- 
pling: a) transitions of the Z^c-system between two multiplets, Rj — > Rj of different dimensionality. 
While these transitions are in general suppressed in the number of colors Nc, a large number of final 
states in higher multiplets can overcome this suppression. The f-channel gluons that induce these 
transitions must be hard to resolve the color structure. Correspondingly, such color excitations can 
be treated perturbatively in an expansion, that is simultaneously a large-A^c/hard scattering expan- 
sion, b) transitions of the Z7c-systems that are rotations within the same color multiplet Ri — > 
Notice that these can involve singlet, as well as octet exchange with a scattering center in the target. 
Such contributions need to be summed to all orders and place no restriction on the hardness of the 
f-channel gluons. 

When viewed from the target side, after the scattering the nucleus will be left in a state with 
multiple color excited nucleons ^ . The higher the number of color excited nucleons, the higher 
will also be the multiplicity of particles produced in the target hemisphere. This clearly implies 
nonperturbative mechanisms for energy loss of the be system. In addition, this picture offers a 
perfect definition of cut Pomerons and connection with the RFT language: each and every color 
excited nucleon must be associated with a unitarity cut of a Pomeron exchanged between beam 



and target. From a technical point of view, the counting of color-octet nucleons is possible after 
a decomposition of the dipole cross section-operator r(^'(C) = r^^^(C) into an 'elastic' 

(4) 

piece that involves two-gluon color singlet exchange with the nucleon, and an excitation piece 
Lex that involves color-octet exchange with a nucleon . We come to the implications on the 
unitarity cut interpretation of the nuclear unintegrated glue. Firstly, one would easily identify the 
cross section for 7-fold quasielastic scattering of a parton a as dOQ^i/d^p = Ca/2CFf^-'\p)- The 
familiar expansion of the nuclear unintegrated glue over multiple convolutions now takes a new in- 
terpretation as an expansion over unitarity cuts corresponding to final states with j color excited nu- 
cleons: 0(b,x,p) =Y,Wj{VAO^))dGQ}i/GQeid^P- In the alternating sign expansion exp[—VA(b)] = 
l)'^/fc![7^/2]'^[/ fi?^p/(p)]'^ of the exponential entering wj, the ^-th order term corresponds pre- 
cisely to an absorptive correction from k uncut Pomeron exchanges. This suggests an interpretation 
of 0(b,x,p) as a coupling of the absorbed gluon-Reggeon to the nucleus. Out of the many possible 
applications we give a brief description on the RFT-cut interpretation of forward-dijets in inelastic 



^Note a certain similarity to tlie concept of wounded nucleons ^3j\. 

^Such decomposition is clearly infrared sensitive, but so is also the transition from cut Pomerons to final state 
particle multiplicities. 
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DIS - the counterpart of the free nucleon result eq.(|l.l[). We obtain the nonlinear -factorization 

for topological cross sections of DIS followed by color excitation of v nucleons (v cut pomerons) 
6. 

da,{tA ^ {qqW TA{h) f ..^ ,2.. ,2 



d^hdzd^pd^-K {InY Jo 



J dp J d^Kid^K2d^Kd{K— Ki — K2 — k) 



dOQd{K) 



d^K 



i/a(j8;z,p- JCi) - VA(j8;z,p- JCi - k) 



2 



X £ 5{v-\-j-k)wj({\-P)v^{h))wMl-l5)v^{h)) Jtj - Kj - (4-1) 

The color channel coupling effects manifest themselves in the presence of two different types of 
cut Pomerons, here they show up as dOQei{K)/d^K oc /(fc), which derives from the singlet-to- 
octet transition of the ^^-pair, whereas dOQj,j{K\) / OQgjd^Ki ■ dOQ^i{K2) /(7Qeid^K2 derive from 
color rotations of the g^-pair in the color-octet state; quark and antiquark rescatter independently 
and contribute j, respectively k color excited nucleons to the final state. The RFT structure is 
further explained in Fig. ^. When going over to topological cross sections for single-jet production, 
say, the single quark spectrum, interactions of the spectator anti-quark do not cancel but leave 
a Cheshire Cat Grin (CCG): the final state will contain color excited nucleons that derive from 
interactions of the antiquark. Amazingly, a novel multiplicity resummation, removes the CCG. 
The contribution from k cut Pomerons can be isolated from 

d^hdzd^p d^hdzd^p 

which roughly corresponds to a re-summation over backward multiplicities hb > k{nB) pp. 

As a consequence of two types of cut Pomerons, common doctrine from hadronic Glauber- 



type models [24] is entirely inapplicable in the pQCD realm. For example an expansion of the 



inelastic Glauber profile function for color dipoles 

r('-') = i_exp[-a(^,r)rA(b)]=exp[-a(^,r)r^(b)]£l[a(^,r)rA(b)]^ (4.3) 

simply has no sensible unitarity interpretation at all. In particular it does not allow one to extract 
topological cross sections from the inclusive inelastic profile function. Additionally, relations on 
the double Pomeron contributions that are in widespread use as a unitarisation snake-oil, such as 
the cancellation of discontinuities 

A2r('") (I^ IP) + 2A2r('"' (1^ y ) = (4.4) 

and the "AGK-ratios" 

A2rD(lPIP) : A2rj"(y IP) : AV^{J^ 1^ ) = 1 : -4 : 2 , (4.5) 



^Here j3 is a dimensionless longitudinal depth of the nucleus, which arises because of the noncommutativity of 
and Sf^', and i/a(j8, ...) denotes a LCWF coherently distorted over a slice of length j8. 
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Figure 2: One of four diagrams representing the RFT unitarity cut for inelastic DIS in terms of two types 
of cut Pomerons. The curly lines represent gluons from the cut Pomeron ly^ associated with the singlet 
to octet transition of the qq pair, and couples to the target through the free-nucleon glue /. Zigzag-lines 
are associated with two cut nuclear Pomerons IF associated with color rotations in the q'^-octet channel. 

/ r,A 

Their coupling through the target is through the nuclear unintegrated glue <I>. When expanded out they give 
rise to the factors °^ dGQ^i{Ki) ■dGQ^i{K2) in eq. (4.1 ). Also indicated are the coherent distortions of the 
Y* qq LCWF through uncut Pomerons. Notice the cylinder/Mandelstam structure of the cut. The coupling 
of multiple Pomerons is local in rapidity and involves no branching tree of triple Pomeron interactions/fan 
diagrams. 



take no account of the intricacies of color channel coupling (two types of cut Pomerons) and the 
subtleties of spectator interactions. See [16] for the -infrared sensitive- relations that replace them. 
Gribov's unitarity relation A2r,o,(IPlP) = —A2roiffractive(^^) between the shadowing correction 
to the total cross section and the diffractive profile function is model independent, and remains of 
course valid. 



5. Summary 

Hard scattering in a regime of strong opacity requires a novel, nonlinear ^x-factorization. Here 
we would not reiterate our points which we made clear enough in the main text, rather we end with 
an outlook: there remains a long shopping list of phenomenological applications connected to the 
production of forward high px-particles/jets, such as: nonperturbative and perturbative quenching, 
its centrality dependence, forward-backward correlations between hard production in the beam- 
and multiplicities in the backward hemispheres, azimuthal decorrelations, etc, and we look for- 
ward to practical applications. On the theoretical side a study of multipomeron couplings and 
the RFT-interpretation of the small-;c evolution of the coherent nuclear glue are to be ordered for 
Diffraction2008. 

It is a pleasure to thank the organisers, especially Igor Ivanov and Alessandro Papa, for their 
efforts, and for invitation to this exciting workshop. 
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